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ABSTRACT 

Context. We consider the Coma cluster of galaxies as a gravitationally bound physical system embedded in the perfectly 
uniform static dark energy background as implied by the ACDM cosmology. 

Aims. We ask if the density of dark energy is high enough to affect the structure of a large rich cluster of galaxies? 
Methods. We use recent observational data on the cluster together with our theory of local dynamical effects of dark 
energy, including the zero-gravity radius Rzg of the local force field as the key parameter. 

Results. 1) Three masses are defined which characterize the structure of a regular cluster: the matter mass Mm, the 
dark-energy effective mass Mde (< 0) and the gravitating mass Mg (= Mm + Mdb). 2) A new matter density profile 
is suggested which reproduces well the observational data for the Coma cluster in the radius range from 1.4 Mpc to 14 
Mpc and takes into account the dark energy background. 3) Using this profile, we calculate upper limits for the total 
size of the Coma cluster, R < Rzg ~ 20 Mpc, and its total matter mass, Mm Ss Mm(Rzg) = 6.2 x 1O 15 M0. 
Conclusions. The dark energy antigravity affects strongly the structure of the Coma cluster at large radii R <; 14 Mpc 
and should be taken into account when its total mass is derived. 
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1. Introduction 

The Coma Cluster of galaxies (A1656) is the most mas- 

■ sive well-studied regular gravitationally bound aggregation 
of matter in the observable Universe. In his classic work 
Zwicky (1933, 1937) applied the virial theorem to the clus- 

■ ter and showed that dark matter dominated the system. 
Zwicky estimated its mass as 3 x 1O 14 M , when normal- 

• ized to the Hubble constant h = 0.71 used hereafter 

Decades later, The & White (1986) found an order of 
magnitude larger value, 2 x 1O 15 M0, with a generalization 
of the virial theorem which must be used when the observed 
sample does not include the entire system. Hughes (1989, 

1998) obtained a similar value (1 - 2) x 10 15 M Q with X-ray 
data under the assumption that the cluster hot intergalac- 
tic gas is in hydrostatic equilibrium. In a similar way, Colles 
(2006) reports the mass 4.4 x 10 14 M & inside the radius of 
1.4 Mpc. A weak-lensing analysis gave 2.6 x 10 15 M Q (Kubo 
et al. 2007) within 4.8 Mpc radius. Geller et al. (1999, 2011) 
extended mass estimates to the outskirts of the cluster us- 
ing the caustic technique (Diafcrio & Geller 1997, Diaferio 

1999) and found the mass 2.4 x 10 15 M Q within the 14 Mpc 
radius. Here the 2a error is 1.2 x 1O 15 M , hence it does not 
contradict the apparently larger mass within 4.8 Mpc. 

In this paper, we re-examine the matter mass estimation 
of the Coma cluster using the data above and a new theory 
model which describes the cluster as a bound spherical sys- 
tem embedded in the cosmic background of dark energy. We 
find that dark energy affects strongly the cluster structure 
at large distances R > 14 Mpc from the cluster center and 
must be taken into account in the matter mass estimate. 



Basic theory is outlined in Sect. 2, three characteristic 
masses of a regular cluster are introduced in Sect. 3, a new 
matter mass profile is defined in Sect. 4, the upper bound on 
the total size of the Coma cluster is calculated and discussed 
in Sect. 5, and the results are summarized in Sect. 6. 



2. Local dynamical effects of dark energy 

We adopt the ACDM cosmology which identifies dark en- 
ergy with Einstein's cosmological constant A and treats 
it as a perfectly uniform vacuum-like fluid with the con- 
stant density pde = 0.71 x 10 -29 g cm~ 3 . The dark energy 
background produces antigravity which is stronger than the 
matter gravity in the present Universe as a whole. This 
makes the cosmological expansion accelerate as discovered 
by Riess et al. (1998) and Perlmutter et al. (1999). 

The cosmic antigravity can be stronger than gravity not 
only globally, but also locally on the scales of ~ 1 — 10 
Mpc (Chernin et al. 2000, 2006; Chernin 2001; Byrd et al. 
2007, 2012), as studied using the HST observations made 
by Karachentsev's team (e.g., Chernin et al. 2010, 2012). 

The local weak-field dynamical effects of dark energy 
can be adequately described in terms of Newtonian mechan- 
ics (e.g., Chernin 2008). Such an approach borrows from 
General Relativity the major result: the effective gravitat- 
ing density of a uniform medium is given by the sum 

p off = p + 3P, (1) 
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where p and P are the fluid's density and pressure (c = 1 
hereafter). In the ACDM model, the dark energy equation 
of state is Pde = — Pde, and its effective gravitating density 



PDEcff = Pde + 3Pde = — 2pde < 



(2) 



is negative, producing antigravity. Einstein's "law of univer- 
sal antigravity" says that a point mass M within uniform 
dark energy generates an acceleration a(r) which includes 
in addition to the Newtonian term a^(r) = —GM/r 2 the 
antigravity effect of dark energy 



geO) = 



4ttG , . 2 8ttG 

— 5— PuEesr /r = H — — 



PDE r - 



(3) 



Then a test particle at the distance R from the center of a 
spherical matter mass Mm (beyond the mass) has the radial 
acceleration in the reference frame of the mass center: 

a(R) = a N (R) + a E (R) = -G^rf + ^-p DE R- (4) 

Eq.4 comes from the Schwarzschild-de Sitter spacetime in 
the weak field approximation. The net acceleration a(R) 
is zero at the distance (Chernin et al. 2000, Chernin 2001, 
2008) 



R - R ZG - [ g7r M 
— Pde 



Gravity dominates at distances R < Rzg , while antigravity 
is stronger than gravity at R > Rzg- A gravitationally 
bound system with the mass Mm can exist only inside its 
zero-gravity sphere of the radius Rzg ■ 

3. Three masses of a regular cluster 

The presence of dark energy in the spherical volume of a 
regular cluster like Coma may be quantified by its effective 
gravitating mass within a given clustrocentric radius R: 

M DE (R) = -yp DE i? 3 = -0.85 x io 1B [^-] 3 M (6) 

The matter (dark matter and baryons) content of the clus- 
ter is characterized by the mass Mm(R) inside radius R: 



M M (R) = An J p(R)R 2 dR. 



Here p(R) is the matter density within the cluster. The sum 

M G (R) = M M + M DE (R) , (8) 

is the total gravitating mass within the radius R. It is this 
mass that can be directly measured by the methods cited in 
Sec.l, which are all related to gravitation and give the grav- 
itating mass Mq(-R), rather than the matter mass Mvr(-R) 
which, however, can be derived from the data of Sec.l using 
Eq.8: M M = M G - M D e- So one has for R = 1.4 Mpc: 

M DE = -2.3 x 10 12 M o , M M ^ M G = 4.4 x 1O 14 M . (9) 

For R = 4.8 Mpc: 

M DE = -9.4 x 10 13 M Q , Mm — Mq = 2.6 x 1O 15 M .(1O) 

And for R = 14 Mpc: 

M DE = -2.3 x 10 15 M o , M M = 4.7 x 10 15 M Q . (11) 



As we see, in the inner cluster at R = 1.4 and 4.8 Mpc 
the dark energy contributes practically nothing compared 
to the gravitating mass, so M G — Mm here. But (curiously) 
the absolute value of the dark energy mass Mde nearly 
equals the gravitating mass M G at R = 14 Mpc; as a result, 
the matter mass Mm — 2Mg — 2|Mde| at this radius. 
The difference between the Eq.ll estimate and the observed 
value of M G is at the level of 4<r (la w \M G , Sec.l). 

In the outer regions of the Come cluster, 

I Mde I > M G , R > 14 Mpc, (12) 
where the antigravity effect thus is significant. 

4. Matter mass profile 

Our estimate of the Coma matter mass within R = 14 Mpc 
(Eq.ll) may be compared with estimates following from 
traditional matter density profiles for dark haloes. 



4.1. NFW and Hernquist profiles 

The popular NFW profile (Navarro et al. 2005) is 

4/5 s 



P = 



(13) 



where R is again the distance from the cluster center, 
p s = p(R s ) and R s are constant parameters. At small radii, 
R « i? s , the matter density goes to infinity, p oc 1/R as R 
goes to zero; at large distances, R>> R s , the density slope 
is p cx 1/R 3 . With this profile, the matter mass profile is 



M M (R) = 167rp s i?3[ln(l + R/R s ) - 



Rj R s 



1 + R/R s 



(14) 



To find the parameters p s and R s , we use the small- 
radii data of Sec.l: M x = 4.4 x 10 14 M Q at R x = 1.4 Mpc, 
M 2 = 2.6 x 10 15 M Q at R 2 = 4.8 Mpc. At these radii, the 
gravitating masses are nearly equal to the matter masses 
there (c.f. Sect.3). The values of M x ,Ri and M 2 ,i? 2 to- 
gether with Eq.14 lead to two logarithmic equations for the 
two parameters of the profile, which can easily be solved: 
R s = 4.7 Mpc, p s = 1.8 x 10~ 28 g/cm 3 . Then we find 
the matter mass within R = 14 Mpc: 



( 7 ) M M ^ 8.7 x 10 15 M c . 



, (15) 

considerably larger (over 70%) than given by Eq.ll. 

Another widely-used density profile (Hernquist 1990) is 



p(R) oc 



R(R + a) 



(16) 



Its small-radius behavior is the same as in the NFW pro- 
file: p — ?* 00 , as R goes to zero. The slope at large radii is 
different: p oc 1/R 4 . The corresponding mass profile is 



Mm(-R) = Mq[ 



R 



R + a' 



(17) 



The parameters M and a can be found from the same data 
as above on Mi, Ri and M 2 , i? 2 : M = 1.4 x 10 16 M Q , a = 
6.4 Mpc, giving another value for the mass within 14 Mpc: 



M M = 6.6 x 10 Mq , R = 14 Mpc, (18) 
Now the difference from the figure of Eq.ll is about 40%. 
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R (Mpc) 

Fig. 1. Matter mass M vs. radius R for three profiles. 1) 
NFW mass profile of Eq.14: M = 13.8[ln(l + R/4.7) - 
4 ^ +R \ - dashed line; 2) Hernquist profile of Eq.17: M 

I ', ^ 6 ^tr ^ - t mn hne; 3) New mass profile of Eq.19: 

M = 8.7 ~ thick hne. Here R is in M P C and 

M is in 10 15 M Q . The values of the matter mass at radii 
R = 1.4,4.8, 14 Mpc are also showed as given by Eqs.9-11 
with la error bars. 




Fig. 2. Matter density p vs. radius R. 1) NFW density pro- 
file (Eq.13): p = R (iJj\ R yi - dashed line; 2) Hernquist pro- 
file (Eq.16): p = wg 9 J + R) 3 ' t nm hne; 3) New density profile 



(Eq.20):p: 



34 

(2.4+K) 4 



thick line. 



4.2. A modified mass profile 

In a search for a more suitable mass profile for the Coma 
cluster, we try a modified version of Hernquist's relation: 



M M (R) = M* 



R 



R + R* 



(19) 



The power of 3 is used now instead of the power of 2 in 
Eq.18. This mass profile comes from the density profile: 



p = —M*RJR- 

47T 



(20) 



The density goes to a constant as R goes to zero; at large 
radii p oc 1/R , like in Hernquist's profile. 

The parameters M* and R* are found again from the 
data for the radii of 1.4 and 4.8 Mpc: M* = 8.7 x 
10 15 M , i?* = 2.4 Mpc. The new profile leads to a smaller 
matter mass at 14 Mpc: 



M 



M 



5.4 x 10 15 M, 



(21) 



which is equal to the Eq.ll value within 15% accuracy. The 
three profiles are graphed in Figs.l and 2. 

5. Discussion 

We now discuss some implications of the above results. 



5.1. Upper limits to size and mass 

The strong effect of dark energy at large radii puts an abso- 
lute upper limit on the total size of the cluster: the system 
can be gravitationally bound only if gravity dominates in 
its volume (as we mentioned in Sec. 2). In terms of the three 
different masses, this criterion may be given in the form 



M G > 0, M M > \M] 



'DE 



(22) 



Both inequalities are met, if the system is not larger than 
its zero-gravity radius ( Eq.5): R < i? max = i?zc- 
If the radius of a system with matter mass Mm is equal 
to the maximal radius R = -R max , its mean matter density 
(see Bisnovatyi-Kogan & Chernin 2012) is 



(pM) 



Mi 



M 



4tt r>3 
3 n ZG 



2/ODE- 



(23) 



This relation and the new profile (Eq.19) now lead to i? max 
and the corresponding matter mass, M max = Mivi(i?m ax ): 

flmax = Rzg = 20 Mpc, M max = M M (i?z G ) = 6.2 x 10 15 M Q . (24) 

For comparison, the other profiles lead to 

i? max = 25Mpc, M max = 1.5 x 10 16 Af Q . (NFW) (25) 

R max = 22 Mpc, M max = 9.1 x 1O 15 M . (Hernquist) (26) 
5.2. Close to the maximal size? 

Our studies of nearby systems like the Local Group and the 
Virgo and Fornax clusters (e.g., Chernin et al. 2010, 2012a) 
suggest their sizes are not far from the zero-gravity radii. 
Around them, flows of galaxies are seen (Karachentsev et al. 
2009; 2010, Nasonova et al. 2011), and the systems are lo- 
cated in the gravity-dominated regions (R < Rzg) and the 
outflows are at R > Rzg- If the local systems have nearly 
maximal sizes, this may explain the apparent underdensity 
of the local universe (Chernin et al. 2012b). 

It is tempting to ask if the matter distribution could ex- 
tend to somewhere near the maximal distance of 20 Mpc in 
the Coma cluster as well. If so, its mass would be near the 
upper limit evaluated above, and still consistent with the 
theory of large-scale structure formation which claims the 
range 2 x 10 15 < M < 10 16 M Q for the most massive bound 
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objects in the Universe (Holz & Perlmutter 2010, Busha 
et al. 2005). Another implication is the predicted (eq. 23) 
mean matter density of the system = twice the dark energy 
density, which does not depend on the density profile as- 
sumed (Merafma et al. 2012; Bisnovatyi-Kogan & Chernin 
2012). Its observational confirmation would directly indi- 
cate the key role of dark energy for the structure of the 
system. Using the cosmological matter density parameter 
Q m = 0.27, its mean density contrast would be 

s= {p)-p m = 2_n^_ 1 = 42 (27) 

Pm * 'm 

Another general prediction (Chernin et al. 2006) is that 
at distances R > Rzq any galaxy m the outflow should 
have a velocity higher than 

U osc = H A R[1 + 2(R ZG /R) 3 - 3{R Z q/R) 2 ] 1/2 (28) 

Here H A = [^pde] 1/2 = 61 kms^Mpc" 1 depends on the 
dark energy density only. Furthermore, within the simpli- 
fied model of Einstein-Straus vacuole, one expects the flow 
to reach the global Hubble rate at the edge of the vacuole 
(« 1.7i? ZG , Teerikorpi & Chernin 2010; Hartwick 2011). 

The situation is complicated by the fact that the Coma 
cluster is not isolated, but lies within the CfA Great Wall. 

5.3. Dark energy estimator 

We have assumed that the dark energy density inside the 
Coma cluster is equal to its global value. One can reverse 
the argument and consider the local dark energy density as 
an unknown constant p x . Its value may be independently 
estimated using our concept of three cluster masses. 

The data (Geller et al. 1999, 2011) give the gravitat- 
ing mass Mq within R — 14 Mpc. The mass Mm at the 
same radius may be found by extrapolating the data from 
R = 1.4 Mpc and R = 4.8 Mpc using a reasonable matter 
density profile (Sect. 4). Then with the masses Mq and Mm 
known for R = 14 Mpc, we may find the dark energy mass, 
Mde(-R) = Mq — Mm, at the same radius. Finally, the dark 
energy density inside the cluster is estimated from 

p x = \M G - M M |g^; i?=14M P c. (29) 

With the density profiles of Eqs. 13 and 16 in Sect. 4, 
Eq.29 gives for the local dark energy density the values 

Px = (1.2 - 2) x 10~ 29 g cm" 3 , (30) 

equal to the global value pde within an order of magnitude 
accuracy (to avoid circular argumentation, we do not use 
the profile of Eq.20 here, as it was partly suggested from 
considerations related to the global value itself). 

6. Conclusions 

1 . Three masses which characterize the structure of a regu- 
lar cluster (like Coma) are defined as functions of the radius 
R: the matter (dark matter and baryons) mass Mm(-R), the 
dark-energy effective gravitating mass A/de (negative) , and 
the total gravitating mass Mq{R) — Mm + Mde- 

2. Of these masses, only the gravitating mass Mq reveals 
itself directly in observations at various distances from the 



cluster center. The dark energy mass Mde may be derived 
using the known global value of the dark energy density. 

3. The mass 2.4 x 1O 15 M measured at R = 14 Mpc by 
Geller et al. (1999, 2011) is the gravitating mass Mq inside 
this radius. The corresponding matter mass is Mm — 2Mq. 

4. At small radii, R << 14 Mpc, dark energy effects 
are almost negligible, |A/de << Mm, and the gravitating 
mass Mq is practically equal to the matter mass Mm- At 
large radii R > 14 Mpc, the antigravity effects are strong 
and | M DE | > Mq. 

5. A new matter mass profile for the Coma cluster re- 
produces well the observational data and accounts for the 
dark energy effects in the radius range from 1.4 to 14 Mpc 
and beyond: Mm(-R) = M*[ r+r ] 3 ; where the constants 
M* and R* can be found from the data for small radii. 

6. The available observational data and the new mass 
profile give upper limits for the Coma cluster total size, 
R <> 20 Mpc, and total matter mass, M M ^ 6.2 x 10 15 M Q . 
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